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ABSTRACT 



We report the results of spectroscopic mapping observations carried out toward 
a ~ I'x 1' region within the northern lobe of the outflow from NGC 2071 using the 
Infrared Spectrograph (IRS) of the Spitzer Space Telescope. These observations 
covered the 5.2—37 fim spectral region and have led to the detection of a num- 
ber of ionic, atomic, and molecular lines, including fine-structure emission of Si"*", 
Fe+, S++, S, the S(0)-S(7) pure rotational lines of H2, the R(3) and R(4) transi- 
tions of HD, and at least 11 transitions of H2O. In addition, the 6.2, 7.4, 7.6, 7.9, 
8.6 and 11.3 jim R\H emission bands were also observed and several transitions 
of OH were tentatively detected. Most of the detected line transitions were strong 
enough to map including, for the first time, three transitions of hot H2O. We find 
that: (1) the water emission is extended; (2) the extended emission is aligned with 
the outflow; and, (3) the spatial distribution of the water emission generally follows 
that observed for H2. The Spitzer/IRS data suggest that both dissociative and non- 
dissociative shocks are present in the portion of the outflow observed here. Most 
of H2, H2O S(2)-S(7), HD, and S I emission detected likely arises from behind the 
non-dissociative shock, while the Si II, Fe II, and most of the H2 S(l) emission orig- 
inate from behind the dissociative shock. Based on the measured line intensities, 
we derive an HD abundance relative to H2 of 1.1 — 1.8x10^ and an H2O num- 
ber density of 12 — 29 cm^'^. The H2 density in the water-emitting region is not 
well constrained by our observations, but is likely between 3x10"^ and 10^ cm~^, 
yielding an H2O abundance relative to H2 of between 2x 10^'^^ and 6x 10^ ^ Future 
observations planned for the Herschel Space Observatory should greatly improve 
the density estimate, and thus our knowledge of the H2O abundance, for the water- 
emitting regions reported here. Finally, we note a possible departure from the H2O 
ortho-to-para ratio of 3:1 expected for water formed in hot post-shocked gas, sug- 
gesting that a significant fraction of the water vapor we detect may arise from H2O 
sputtered from cold dust grains. 

Subject headings: ISM: Abundances, ISM: Jets and Outflows, ISM: Molecules, 
Stars: Formation, Stars: Winds, Outflows 
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1. Introduction 



Outflows from young steUar objects provide a natural laboratory toward which it is pos- 
sible to test the predictions of chemical models thought to apply to warm (T > 300K), dense 
(n(H2) ^ 10^ cm"'^) gas. Among these predictions is that water will be for r ned rapidly in the 



warrn gas via a set of neutral-neutral reactions (cf. Elitzur & de Jong 



1978 



Elitzur & Watson 



19781). Moreover, these reactions can potenti ally convert all of the oxygen not locked in CO into 
H2O (cf. iBergin. Melnick. & Neufeldlll998h . Indeed, the gas-phase water abundances toward 
shock-heated gas ass ociated with outflows is found to be higher than that within cold, quiescent 
molecular clouds (cf. 



20001 : 



Harwit et al. 



Benedettini et al. 



1998 



Melnick et al. 



2000; 



Neufeld et al.l 



20001 : 



Nisini et al. 



2002j), although the inferred abundances rarely approach [H20]/[H2] 



10~^ as would be expected for the full conversion of all oxygen not in CO into H2O. In addition, 
an unknown fraction of the elevated gas-phase water abundance is likely due to the removal of 
water -ice from grains that pass through nondissociative shocks (|Draine. Roberge. & Dalgamo 
19831). Based on a detailed analysis of ground-state ortho-water spectra obtained by the Sub- 
millimeter Wave Astronomy Satellite (SWAS) toward 18 outflow sources. Franklin et a/.l ( 2008 ) 
propose that the lower-than-expected water abundances within these sources is due to the small 
fraction of outflow material that actually passes through strong shocks. This also implies that 
the beam filling factor of regions with highly-abundant water is lower than previously assumed. 

Unfortunately, these previous observations lacked the spatial resolution to determine the 
detailed water distribution. Thus, the interpretation was hampered by uncertainties in both the 
fraction of the beam filled with H2O emission and the physical conditions in the H20-emitting 
gas as might be revealed by co-spatial emission from H2. With beam sizes that range between 
1.8" and 5.5", depending on wavelength, and high sensitivity, the Infrared Spectrometer (IRS) 
aboard the Spitzer Space Telescope has enabled the study of H2O, H2, and HD with unprece- 
dented detail. In this paper, we describe observations of the outflow source NGC 2071 and 
report the first detection of spatially extended, hot H2O vapor with Spitzer/IRS. In addition to 
H2O, we have also mapped the ground-vibrational state distribution of H2, as well as low-lying 
rotational and fine-structure emission from HD, S I, Fe II, and Si II. 

NGC 2071 is a region of active star formation, located inside the northern part of th e 
Orion B (L 1630) molecular cloud complex at a distance of 390 pc (|Anthony-Twaroglll982|) . 
The region of interest in this study lies approximately 4' north of the optical reflection nebula 
NGC 2071 and is the site of a prominent outfl ow. Fig.JT] sho ws a portion of the outflow, as 
delineated by the H2 1-0 S(l) 2.12/im emission teisloffellboool) . along with the region observed 
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here. Near-infrared continuum imaging and spectroscopy of this ar ea reveal a cluster of infrared 
sources near the center of this outflow dWalther et a/.lll99lL Il993 l). The detection of a bipolar 
thermal radio jet from source IRS3 oriented in the direction of the NGC 2071 N orth flow, and 
of a ring of H2O maser spots perpendicular to the flow axis ( Torrelles et a/.l 1998|) . suggest that 



IRS 3 is the driving source of the NGC 2 071 North fl ow. 
ha s bee n mapped previou sly in CO JsneUg? a/.lll"984l : 



le outflow, which extends more than 



19891), CS (IZhou et q/.II1991I). SO. SiO. and HCO+ (e.g.. 



Moriarty-Schieven. Hughes & Snell 



Chemin & Mas son 



in H2 1-0 S(l) (e.g.. lEisloffel 
as high as 60 km s~^ 



1991 



19931), and 



2OOOI and references therein). The presence of high-velocity gas 



Chernin & Masson 



and the enhanced abundance of species such as SO and SiO (e.g. 



19931) suggest that shock processing of the outflow gas is highly likely. 



In this paper, we present the results of Spitzer/IRS observations toward the northern lobe 
of the NGC 2071 outflow and discuss the new understanding revealed through high spatial res- 
olution mapping in the lines of H2 S(0)-S(7), HD R(3) and R(4), Sill, Fell, SI, and three 
transitions of Hg^^O: 725— 616, 634—505, and 734—625. In addition, we present spectra that in- 
clude several other H2O features along with several lines that we tentatively ascribe to OH. 
Section 2 summarizes the observations and Section 3 presents the results. Section 4 describes 
the calculations used to analyze the line data and the conclusions are discussed in Section 5. 



2. Observations and Data Reduction 

NGC 2071 was observed on 20 and 21 October 2004 as part of Spitzer Space Telescope 
Cycle 1 General Observer (GO) program 3423. The observations utilized the Long-High (LH), 
Short-High (SH), and Short-Low (SL) modules to achieve the highest possible spectral resolv- 
ing power and complete spectral coverage available to IRS. The mapped areas were approx- 
imately 1' X 1' for LH and SH and approximately 1' x 1.5' for SL, obtained by stepping in 
spacings of half- slit width perpendicular to the slit and in ~ 4/5 of the slit length parallel to 
the slit. The raw data were processed with the Spitzer Science Center (SSC) pipeline version 
S12.0.2, which provided basic calibrated data (BCD). Further data reduction utilized a data 
pipeline developed by our group. This pipeline in corporates procedure s for removing bad or 



"hot" pixels, utilizes the SMART software package (|Higdon et a/.ll2004|) for calibration and ex- 
traction of spectra at each position along the slit, applies a slit loss correction function (SLCF) 
to correct the flux calibration performed for point-like sources to extended sources, and auto- 
matically fits Gaussian lines and second-order polynomials to a set of predetermined spectral 
lines at each position. The first three steps of this data processing are described in more de- 
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tail in (INeufeld et a/.ll2006a|) . The line fitting was performed iteratively to achieve the best-fit 
parameters. In some cases, the resulting line maps exhibited noticeable striping attributable to 
varying pixel sensitivity across the slit and by the after-effects of correcting for bad pixels. A 
pixel response correction factor (PRCF) was determined to remove striping in the final maps of 
all LH and SH observations. We used a statistical approach to determine the relative sensitivity 
in each pixel w hile preserving the to tal flux in all pixels across the slit. Details of this method 
are provided in 



Neufeld et al. 



( 2007 ). 



3. Results 

The observed spectra show emission from molecular, atomic and ionic species as well 
as a number of PAH features. In particular, we detect emission from [S I], [S III], [Fell] and 
[Sill] transitions, H2 rotational lines from S(0) to S(7), HD R(3) and R(4) transitions, and at 
least 1 1 H2O transitions. Most of the line emission was strong enough to detect over extended 
regions and these maps are shown in Figs. [2]-|6l All maps have been normalized to the peak 
line intensities listed in Table 1 . For reference, the positions of the three strongest peaks in the 
H2 S(l) 17 jjm emission are marked as PI, P2, and P3 on these maps. As can be seen in Fig.lH 
the strongest water emission lies close to the peak of the H2 S(l) 17 fim emission (position PI 
in Fig. 131). It is also clear that: (1) the water emission is extended; (2) this extended emission is 
generally aligned with the outflow; and, (3) the water emission exhibits secondary peaks toward 
positions of enhanced H2 emission. Fig. |7] shows a direct comparison of the spatial distribution 
of the H2O (725 — 616) emission versus that due to the ro-vibrational H2 1-0 S(l) and pure 
rotational H2 S(l) emission. While the overall distribution of the two species is quite similar, 
the peaks of the H2O and H2 are sometimes offset with respect to each other. For example, 
the peaks of the H2O 725 — 616 and the H2 S(l) emission are offset by 3.8", which exceeds the 
stated IRS sky registration error, and thus is likely real. 

The spectra obtained toward the water emission peaks closest to positions PI and P2 are 
shown in Figs. [8]— [lOl In order to reduce the baseline noise, we averaged the data toward 
the water peaks using a Gaussian weighing function with a wavelength-dependent HPBW of 
between 1 1.0" at 27.0 /im and 14.7" at 5.5 fim. The final synthetic bea ms have effective beam 



widths of 15" for all wavelengths of interest (cf. iNeufeld et a/.ll2006a ). The spectrum labeled 
H2O Peak 1 is centered at RA 05^ 47°^ 04?6, Dec -i-00° 22' 52.2" (J2000), while the spectrum 
labeled H2O Peak 2 is centered at RA 05^ 46" 40?4, Dec -i-00° 22' 43.0" (J2000). For almost aU 
of the H2O and H2 transitions detected, the 15" synthetic beam is sufficiently large to encompass 
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the water emission as well as the bulk of the H2 emission toward each peak. The measured line 
fluxes for the water features associated with H2O Peak 1 are provided in Table 2. Unfortunately, 
at the spectral resolving power of the Spitzer/IRS, many of the H2O features are blends of two 
or more H2O transitions. The extinction-corrected line fluxes for the other key species in this 
spectrum are given in Table 3. We also note the tentative detection of emission from a number 
of OH rotational transitions (see Figs. [8] and |9l) whose upper levels lie as much as ~ 3500 K 
above the ground. 

Direct measures of the extinction toward the NGC 2071 outflow are not available. In- 
stead, we inferred the required extinction corrections as a function of wavelength in the Spitzer 
bandpass by fitting the wings of the silicate absorption feature at 9.7 /im, hence minimizing 
the effects of line saturation when present. Our silicate fitting model uses the renormalized syn- 
thetic Gal actic extinction curve per unit hydrogen column density calculated for Ry = 3. 1 or for 
Ry = 5.5 JDrainell2003albl) in linear combination with a first-degree polynomial and 4 Gaussians 
with fixed widths and fixed centers to account for emission from interfering lines. Since the gas 
clumps probed by our data have dust properties most probably intermediate between diffuse 
and dense molecular clouds, the extinction corrections we adopted represent the mean of the 
extinction corrections we obtained for _Rv = 3.1 and those we obtained for Ry = 5.5. 



4. Water Calculations 



To assess how the H2O line strengths measured here constrain the water abundance, the 
equilibrium level populations of all ortho and para rotational levels of the ground H2^^0 vibra- 
tional state with ene rgies EA up to 7700 K hav e been calculated using an escape probability 
method described by lNeufeld & Melnickl(|l99l[) . In modeling the excitation of water vapor, we 
are limited by the availability of molecular data. Thus far, at the high temperatures of relevance 
here, a com plete set of rate coefficients has only been computed for coUisional excitation of 
H2O by He (iGreen. Maluendes & McLeanlll993b . A new and accurate nine-dim ensional po- 
tential energy surface has recently be en computed for H2O-H2 (iFaure et q/.II2005|) . and quantal 



calculations are currently underway (IDubernet 



2007 



personal communication), to obtain rate 
coefficients for the coUisional excitation of H2O by H2. However, quantal results for transi- 
tions involving the highly-lying rotational states observed by Spitzer have not yet been com- 
pleted. Neverthele ss, for the subset of those transitions that show the largest rate coefficients. 



Faure et al. 



(|2007() have recently performed quasi-classical trajectory (QCT) calculations to es- 



timate the coUisional rates, and have kindly made the results available to us prior to publication. 
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Unfortunately, transitions with small rate coefficients still play a significant role in the exci- 
tation of the s pectral lines that we have observed with Spitzer; for those transitions, we have 
had to use the lGreen. Maluendes & McLeanI (| 1993b results, scaled by a factor 1 .348 to account 



for the difference between the H20-He and H2O-H2 reduced masses. Finally, a polynomial fit 
to the exact expression for the photon escape probability from a plane-parallel emitting region 
(IHummer & Rybickilll982|) is used. 



We assume that the physical conditions within H2O Peak 1 are those required to repro- 
duce the measured H2 S(0) through S(7) flux from the same region. Table 4 summarizes these 
physical conditions. The model H2 column densities are assumed to be those consistent with a 
simple shock model in which 



N{H2) 



6.45 X 10 



20 





0.5 


J- gas 


-0.555 


105 




1000 





cm 



(1) 



Neufeld et al. 



2006a ). We assume that the 



where n{H2) is the density in the postshock gas (cf. 
fraction of the beam filled by the warm and hot H2O is given by the ratio of the measured values 
of N(H2) given in Table 4 to those derived above. To estimate the physical parameters from 
the eight H2 rotational lines, we have performed calculations in which the emitting region is 
approximated by a two-component model for two different cases: (1) clouds with homogenous 
density; and, (2) clouds with homogenous pressure. The six best-fit parameters (temperature, 
column density and H2 ortho-to-para ratio for the two components) are determined by mini- 
mizing the root mean square of In (measured flux/computed flux) for the H2 S(0) to S(7) lines. 
Results were obtained for a series of H2 number densities ranging from 10^ to 10"^ cm^'^, and 
pressures (x \/k) ranging from 10^ to 10^° cm^'^ K. The greatest differences between the two 
models - homogenous density or pressure - appear at low densities (or low pressures), where the 
pressure-balanced assumption requires a relatively higher temperature for the hot component. 
The results of the two models converge as LTE is approached for the higher-lying transitions, 
i.e., for number densities greater than 10^ cm~'^ (or pressures/fc greater than 10^ cm~'^ K). In 
this case, the best-fit temperature for the warm component is around 300-360 K while that for 
the hot component is around 1000-1500 K. 

Using the physical conditions that best fit the H2 observations, the abundances of H2O and 
HD toward H2O Peak 1 were then derived. For H2O, the total (ortho-i-para) abundance at each 
assumed density was varied to produce the best overall fit to the features listed in Table 2. We 
choose these features to fit because they are the strongest observed transitions between 5 and 
36.5 fim toward H2O Peak 1; both the observed and predicted line fluxes for H2O transitions 
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lying at A < 29 fim are less than 5 MJy/sr, and most typically less than about 3 MJy/sr, render- 
ing their unambiguous identification above the noise difficult. As discussed further below, the 
H2O ortho-to-para ratio was treated as a free parameter. Similarly, the abundance of HD was 
determined by fitting the measured R(3) and R(4) line fluxes, although the best-fit HD abun- 
dance is not a sensitive function of number density due to the moderately low critical densities 
of the R(3) and R(4) transitions of ~ 10^ to 10^ cm~'^. Finally, it would be useful to also fit the 
measured OH rotational lines. Unfortunately, the lack of reliable coUisional rate coefficients for 
these high-lying transitions makes the results of such a calculation uncertain at this time. 

As shown in Fig. [TTl fits to the H2 and HD emission lines alone favor gas with higher 
densities or pressures. The best-fit HD abundance relative to H2 is between 1.1x10"^ and 
1.8x10"^, with only a weak dependence on density as noted above. This abundance is very 
similar to those derived previously from o bservations of HD R (3) and R(4) in other shocked 
regions that we have observed with Spitz.er dNeufeld o/l 2006b ). Fits to the H2O features were 
computed in two ways. First, fits to ten of the eleven H2O features in Table 2 were calculated 
assuming an H2O ortho-to-para ratio of 3:1 as well as an ortho-to-para ratio that was allowed 
to vary in order to achieve the best overall fit at each density between 10'^ and 10^ cm The 
strong feature at 34.55 fim was not included. The conditions needed to reproduce the flux in 
this line resulted in fluxes in the other H2O transitions that are in significant disagreement with 
the observations. Thus, we conclude that either this line is blended with an unknown feature 
or uncertainties in the collision rates may be resulting in lower-than-actual populations in the 
upper 734 state. Second, to better isolate the effects of varying the H2O ortho-to-para ratio, 
we examined the three strongest features in which the only contribution to the line flux was 
an ortho-H20 transition, at 29.84 /im, or a para-H20 transition, at 33.13 and 36.21 fim. These 
results are also shown in Fig.[TTJ Based on fits to the water data alone, H2 densities in the water- 
emitting region of between about 3x10^ and 3x10^ cm^'^ appear favored. The best fits to the 
isolated ortho and para features are achieved for H2 densities between 3x10"^ and 10^ cm^^ and 
an assumed H2O ortho-to-para ratio of between 1.5 and 1.6. 

The best-fit water abundances are shown in Fig. [121 To a good approximation, the best fits 
to the H2O line fluxes are proportional to the H2 density x the water abundance, resulting in a 
derived number density of H2O molecules of 12 — 29 cm^^, independent of H2 density. This is 
a consequence of the fact that the observed transitions have critical densities > 10^ cm and 
are optically thin under the conditions considered here. In this limit, the flux in an ortho-H20 
(or para-H20) line is 



Page 9 



where A^e is the transition spontaneous emission rate, is the transition frequency, ri(H2) is 
the H2 density, £ is the depth of the emitting region, x(H20) is the water abundance relative to 
H2, /„ is the fractional population in the upper state of the emitting transition, is the fraction 
of all H2O molecules in the ortho (or para) state, and is the solid angle of the water emitting 
region. 

Fig. [T3] shows the resulting 29 — 36.5 fim water spectrum computed for an H2 density of 
10^ cm~^ and H2O ortho-to-para ratios of 3:1 and 1.55:1, the latter being the best-fit ortho-to- 
para ratio at a density of 10^ cm^'^. The computed spectra were convolved with a Gaussian 
profile whose width was selected to match the IRS Long-High resolving power. A/ A A, of 600 
(FWHM). As can be seen in Fig. [131 the differences between the spectrum computed assuming 
an H2O ortho-to-para ratio of 3: 1 and 1 .55: 1 are small. Nevertheless, assuming an ortho-to-para 
ratio less than 3:1 results in a slightly better overall fit to the H2O data and suggests a possible 
departure from the usual assumption that water in post-shocked gas is always present in the LTE 
ortho-to-para ratio of 3: 1 . This tentative finding is discussed further in the next section. 

Taking into account the H2, HD, and H2O spectra, the overall best fits to these data favor an 
H2 density in the range 3x 10^ to 10*^ cm"^, and 2x 10"^ < x(H20) < 6x 10"^. Unfortunately, 
the present data do not allow us to better constrain the H2 density, and thus the water abundance, 
within the emitting region. However, as we discuss in the next section, future observations 
planned for the Herschel Space Observatory may permit more direct measures of the H2 density. 



5. Discussion 



The capabilities of the Spitzer/IRS instrument has enabled the study of outflow regions 
in ways not possible previously. In particular, the ability to map the distribution and intensity 
of key species, like H2, HD, Fe+, Si+, S, and, of course, H2O, with a spatial resolution of 
a few arcseconds permits more stringent tests of our understanding of the structure and com- 
position of outflows. Earlier space missions, such as the Infrared Space Observatory (ISO), 
SWAS, and Odin have establishe d the correlation be t ween outflows from young stellar objects 
and strong water emission (e.g.. 



Harwit et al. 



1998; Melnick et al. 



2000; 



Nisini et al.llOOd. loiannini. Nisini. & Lorenzettill2001l ; IBenedettini et all 



Neufeld et al. 



2000 : 



20021). Unfortunately, 



these observations were conducted with large beams relative to Spitzer/IRS or, in the case of 
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the ISO Short- Wavelength Spectrometer with its 10" beam, lower sensitivity than Spitzer/IRS. 
As a result, the filling factor of H2O within the beam was not always directly known or key 
species helpful in the interpretation of the water data, like H2 and HD, were not detected. The 
observations reported here remedy most, but not all, of these previous shortcomings. 

The primary goal of this study is the further understanding of water in outflows. Unlike 
the production of water in quiescent molecular clouds, which proceeds via a set of relatively 
slow ion-neutral reactions in the gas phase (e.g. jHerbst. E.. & Klempereilll973[). or through th e 



HoUenbach et al. 



20081), 



photodesorption of water formed on the surface of dust grains (e.g., 
the material in outflows often passes through a shock which both comp resses and he ats the 
gas. In the post-shocked region of non-dissociative, or C-type, shocks (lDrainelll980|) . with 
shock velocities, Vs, ^ 10 km s~^, temperatures exceed 300 K enabling a pair of neutral-neutral 
reactions (H2 + ^ OH -1- H and Hp + OH ^ H^^Q -i- H) tha t can rapidly convert all gas-phase 
oxygen not bound in CO into H2O (lElitzur & de Jonglll978|) . For Vg ^ 50 km s^, most C-type 
shocks break down, giving rise to dissociative, or J-type, shocks. Molecules are completely 
dissociated, either by the precursor UV field generated by the hot gas near the shock front or by 
colli sions in the shock, and then r eform on grain surfaces in the cooler (~400 K) downstream 
gas (|Hollenbach & McKeelll989|) . Because C-type shocks have very low ionic abundances, 
strong Si II (34.8 /xm) and Fe II (26.0 /xm) emission is considered a tracer of J-type shocks. 

The Spitzer/IRS data suggest that both J- and C-type shocks are present in the portion of 
the outflow observed here. The strong Si II and Fe II emission evident in Fig. [6] follows the out- 
flow, peaking close to H2 positions P3 and P2, and lies closer to the source of the outflow than 
the most prominent H2 position, PI. This may indicate the presence of fast-moving gas along 
the center of the outflow which has been slowed as it approaches position P I. This picture is sup- 
ported by previous CO J = 3— 2 observations (|Chemin & Massonlll992|) . which revealed well 
coUimated, high velocity (30 — 60 km s~^) gas in a symmetric bipolar outflow only within 1.5' 
of the source, presumed close to IRS 3. Moreover, the CO data show velocity peaks appearing at 
± 1' from the source. Further support is provided by the pronounced pe ak in the S I (25.2 /xm) 



HoUenbach & McKee 



emissi on (see Fig. [6]) near position PI. In the J-type shock models of 
(119891) . for densities between 10^ and 10^ cm~^, the strongest S I emission arises at the lowest 
velocities (i.e., Vs ~ 30 km s^^), whereas both the Si II and Fe II emission incre ase with shock 



veloc ity. In addition, previous Spitzer/JRS observations of supemovae remnants (|Neufeld et al. 



2007|) show a strong correlation between the SI and H2 S(l)— S(7) spatial distribution, sug- 



gesting that most of the S I emission arises in the slower, nondissociative shocks. This could 
account for the difference in the spatial distribution of Si II, Fe II, and S I. 
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That the H2 S(l) peak emission appears to be located closer to the source of the outflow 
(near IRS 3) than all of the higher-lying H2 transitions, i.e. S(2) through S(7), as well as that of 
the H2O may also be explained by this scenario. In the multiple shock picture, the H2 S(l) emis- 
sion could arise primarily in the ~ 400 K gas behin d a moderate J-type shock - the post shock 
region in which it is predicted that H2 is reformed (|Hollenbach & McKeelll989t) . Meanwhile, 
the higher- lying H2 and H2O lines would arise from behind the C-type shock slightly farther 
from the outflow source in which gas temperatures can heat molecules to more than 2000 K. 

The good fit to the H2 and HD data (see Fig. [TTI) obtained with conditions characteristic 
of a C-type shock provide the best evidence for the coexistence of a non-dissociative shock. 
Based on the morphology of the Sill, Fell, and H2 emission, it is tempting to suggest that 
the C-type shock emission arises in the entrained slower moving shocked material surrounding 
the high-velocity J-type shock gas which, itself, slows and transitions to a C-type shock near 
position PI. 

That the water emission follows closely the H2 emission is in accord with the predic- 
tions of C-type shock models. Unfortunately, our inability to unambiguously determine the H2 
density in the water-emitting region from the Spitzer/IRS data prevents us from fully testing the 
other key prediction of the C-type shock models - namely that the water abundance will be high 
(i.e., x(H20) > 10^^) in the post -shocked gas. Oth er observations provide general support for 



the density range preferred here. IZhou et all (|199l|) have observed the NGC 207 1 outflow in the 



CS J = 2-1, 5-4, 6-5, and 7-6 transitions. Their best fit to the CS wing emission yield densities 
of 1 — 4 X lO^cm^'^. Unfortunately, the varying beam sizes with which these observations were 
conducted, ranging between 1 1" and 24", along with uncertainties in the gas temperature and 
the surface filling factors for the different transitions contribute to the overall uncertainty in the 
density determination. Nonetheless, if this range of densities applies to H2O Peak 1 region, then 
the best-fit water abundance would lie between approximately 5 x 10^^ and 2 x 10~^ relative to 
H2 which still represe nts a significant a bundance enhancement above that found in quiescent 



molecular clouds (e.g. JSnell et a/.ll2000|) and is consistent with the predictions of C-type shock 
models. 

Future observations of far-IR high- J CO lines - planned for the Photodetector Array Cam- 
era and Spectrometer (PACS) instrument aboard the Herschel Space Observatory - promise to 
yield a valuable means of estimating the H2 density in NGC 207 1 . Given the gas temperatures 
and column densities derived from our Spitzer observations of the H2 S(0)— S(7) rotational tran- 
sitions, we have computed the CO line intensities expected for a CO abundance of 10 relative 
to H2. The sample results shown in Fig. [14] were obtained by means of a statistical equilibrium 
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calculation in which we adopted the molecular data described in iNeufeld & KaufmanI (| 19931) . 
Fig. [14] indicates that far-IR transitions of CO should be readily detectable with Herschel and 
that their strengths will provide a key probe of the gas density over exactly the desired range: 
4 < logio(H2/cm^'^) < 6.5. Here, the dashed line corresponds to the Herschel/FACS spectral line 
sensitivity of ~ 2x 10~^ erg cm^^ s^^ sr {5a in 1 hour, corresponding to a flux of 5 x 10^^^ 
W into a 9.7 x 9.7" pixel). Because the CO lines are optically-thin, the predicted line inten- 
sities are proportional to the assumed CO abundance. However, the far-IR line ratios (Fig.[T5l) 
are independent of the assumed CO abundances and would provide a good estimate of densities 
in the range 5 < logio(H2/cm^'^) < 6.5 without any required knowledge of the CO abundance. 

The question of the distribution and abundance of the gas-phase water is important not 
only as the answers test the validity of long-standing chemical predictions. Many models of the 
larger ambient clouds take as a starting point certain molecular enrichment due to th e passage 
of nearby shocks . A recent analysis of 1 8 molecular outflows by Franklin et al\ dlOOSh based on 
observations of the ground-state ortho-H20 transition from SWAS and the J = 1—0 transition 
of ^^CO and ^^CO obtained with the Five College Radio Astronomy Observatory shows that the 
ortho-H20 abundance (relative to H2) in most outflows is between 10^'' and 10"^, assuming the 
H2O and CO emission arises in the same gas. However, an examination of the water abundance 
as a function of outflow velocity reveals a strong dependence; the water abundance rises with 
outflow velocity, reaching abundances of ~ 10^'^ at the highest velocities. However, the mass 
associated with the highest velocity emission is found to be small compared with the total out- 
flow mass, leading to the conclusion that only a small fraction of the outflowing molecular gas 
has passed through shocks strong enough to fully convert the gas-phase oxygen to water. Using 
the measurements obtained here, combined with the data soon to be available from Herschel, it 
will be possible to independently test this conclusion. 

Finally, our potential detection of an H2O ortho-to-para ratio that is not in equilibrium at 
the gas temperature may hint at the origin of the water vapor emission we detect. The lowest 
energy level of para-H20 is ~ 34 K below that of ortho-H20 and, as such, the ratio of the ortho 
and para populations in LTE is temperature dependent. When water forms in the gas phase via 
exothermic reactions the energy released is much greater than this energy difference and the 
ortho-to-para ratio will reflect the 3:1 ratio of statistical weights between these species. When 
water forms on the surfaces of cold dust grains, it is believed that the energy generated by the 



chemical reaction is shared with the grain and the water molecules shou 
ortho-to-para ratio that reflects the grain temperature (e.g.. 



Limbach et al. 



d equ ilibrate to an 



20061) . If the grain 



temperature is below ~ 50 K, then the ortho-to-para ratio will lie below 3:1. In this regard, it has 
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been known for many years that water in numerous cometary comae exhibits an ortho-to-para 
ratio below 3:1, which has been attributed to forma tion at temperatures of ~25— 30 K (e.g., 



Mumma. Weaver & Larson 



1987 



Bonev et al. 



2007 ). 



Our measurements, while uncertain, potentially limit the ortho-to-para ratio to between 
1.5-2.5, or equivalent temperatures of ~ 20—30 K. If the water is produced entirely in shocks 
via the well known sequence of neutral-neutral reactions, the ortho-to-para ratio should be 3: 1 . 
However, if the water enters the gas phase via sputtering from cold grain surfaces, then the ratio 
would reflect equilibrium at the grain temperature. A mixture of both sputtering and gas-phase 
formation cannot be ruled out and thus this measurement represents an upper limit to the ortho- 
to-para ratio of water ice. If our measurement is confirmed, this would therefore suggest that a 
significant amount of the H2O emission we detect has its source in water vapor produced from 
grain sputtering and not via gas-phase reactions. 
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Table 1. Peak Intensity for Each Mapped Transition 






Rest 


Peak 


Species 


Transition 


Wavelength 


Intensity 






/im 


10 ^ erg cm ^ s ^ sr ^ 


H2O 


725 - 616 


29.837 


4.39 




634 - 5o5 


30.899 


7.08 




854 - 743 


30.871 






734 - 625 


34.549 


7.08 


H2 


8(0) 


28.221 


21.10 




S(l) 


17.035 


189.31 




S(2) 


12.279 


1721 




S(3) 


9.665 


1421 






8.025 


1200 




S(5) 


6.909 


3300 




S(6) 


6.109 


1746 




8(7) 


5.511 


1800 


HD 


R(3) 


28.502 


3.37 




R(4) 


23.034 


6.20 


SI 


'Pi - 'P2 


25.249 


201.70 


Fell 


0^07/2 - a^Dg/s 


25.988 


40.84 


Sin 


^P3/2 - 


34.815 


55.41 
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Table 2. 


Beam-Averaged Line Fluxes for the H2O Features Toward H2O Peak 1 










Extinction-Corrected 


TT 

H2O 








T inp FliiY V 1 0"^-'^ 


Feature 


Para (p) 


Transition 


Wavelength 


into 15" Beam^ 










Wcm-2 


36.21 


P 


624 - 5 15 


36.212 


1.57 


35.92 





652 - 541 


35.938 


0.70 




P 


651 - 542 


35.904 




35.67 





845 - 7,34 


35.669 


0.55 


J J Ad 


P 


J33 - '^■04 


c /no 

35.472 


2.43 







'43 - O34 


35.428 




34.55 





^34 - 625 


34.549 


2.02 


33.13 


P 


753 - 642 


33.127 


0.64 







I61 16 — 15o 15 


33.042 


1.42 




P 


I60I6 ~ 15ll5 


33.042 









661 - 550 


33.005 






P 


660 - 551 


33.005 









'52 - O43 


32.991 









^-'114 ^^213 


32.965 






P 


15214 -14i 13 


32.960 




31.75 





441 - 3i2 


31.772 


0.70 




P 


844 - 735 


31.738 




30.89 





634 - 5o5 


30.899 


1.82 







854 - 743 


30.871 




30.53 


P 


762 - 651 


30.529 


0.63 







^61 ~ 652 


30.526 




29.84 





725 - 616 


29.837 


1.27 



^ At the spectral resolving power of Spitzer/IRS, some H2O features may 
be blends of two or more H2O transitions. The line fluxes listed here rep- 
resent the total flux for each H2O feature. 
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Table 3. Beam- Averaged Line Fluxes for Key Species Toward H2O Peak 1 



Extinction-Corrected 
Rest Line Hux x 10"^^ 



Species 


Transition 


Wavelength 


into 15" Beam 








w cm 


H2 


S(7) 


5.5112 


275.2 


H2 


S(6) 


6.1086 


110.9 


H2 


S(5) 


6.9095 


513.0 


H2 


S(4) 


8.0251 


207.0 


H2 


S(3) 


9.6649 


474.4 


H2 


S(2) 


12.2786 


138.1 


H2 


S(l) 


17.0348 


68.3 


HD 


R(4) 


23.0338 


2.2 


OH 


2n3/2 23/2+ - 21/2- 


25.0351 


0.41 


OH 


^Uy2 23/2- - 21/2+ 


25.0899 


0.48 


[SI] 


'Pi - ='P2 


25.2490 


78.8 


[Fell] 


0,607/2-0,609/2 


25.9883 


7.6 


OH 


2n3/2 21/2- - 19/2+ 


27.3935 


0.13 


OH 


21/2+ - 19/2- 


27.4546 


0.58 


OH 


Tli/2 19/2 - 17/2+ 


27.6516 


0.31 


OH 


2ni/2 19/2+ - 17/2- 


27.6967 


0.26 


H2 


S(0) 


28.2188 


9.1 


HD 


R(3) 


28.5020 


2.3 


OH 


2n3/2 19/2+-17/2- 


30.2772 


0.25 


OH 


2n3/2 19/2--17/2+ 


30.3459 


0.60 


OH 


2ni/2 17/2+-15/2- 


30.6573 


0.33 


OH 


2ni/2 17/2--15/2+ 


30.7063 


0.37 


[sm] 


'Pi - 'Po 


33.4800 


1.5 


[Sill] 


^P3/2 - ^Pl/2 


34.8152 


18.7 
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Table 4. Best-Fit Physical Conditions Within HgO Peak 1 Inferred from H2 S(0)-S(1) Lines 



Logio H2 Density (cm ^) 



Parameter 


4.0 


4.5 


5.0 


5.5 


6.0 


6.5 


7.0 


Warm Component: 
















Temperature (K) 


377 


361 


352 


344 


336 


330 


328 


Logio N(H2) (cm~^) - measured . . . 


20.60 


20.61 


20.61 


20.60 


20.59 


20.59 


20.58 


Logio N(H2) - theory 


20.54 


20.81 


21.06 


21.32 


21.57 


21.83 


22.08 


H2 Ortho-to-Para Ratio 


0.69 


0.71 


0.68 


0.63 


0.58 


0.55 


0.53 


Beam Filling Factor 


1.0 


0.64 


0.35 


0.19 


0.10 


0.06 


0.03 


Mass of Warm H2 (x 0.1 M©) 


0.04 


0.07 


0.12 


0.21 


0.38 


0.68 


1.22 


Hot Component: 
















Temperature (K) 


2111 


1621 


1341 


1169 


1070 


1024 


1006 


Logio N(H2) (cm~^) - measured . . . 


19.31 


19.38 


19.47 


19.58 


19.66 


19.71 


19.73 


Logio N(H2) (cm~^) - theory 


20.13 


20.44 


20.74 


21.02 


21.29 


21.55 


21.81 


H2 Ortho-to-Para Ratio ^ 


3.45 


3.58 


3.68 


3.68 


3.62 


3.56 


3.53 


Beam Filling Factor 


0.15 


0.09 


0.05 


0.04 


0.02 


0.01 


0.008 


Mass of Hot Ho (x 0.1 M;) 


0.01 


0.03 


0.06 


0.11 


0.20 


0.36 


0.65 



Errors on the H2 ortho-to-para ratio are hard to estimate, but the best-fit values for the hot 
component are probably consistent with the LTE value of 3. 
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Fig. 1. Map of the northern lobe of the ou tflow from NGC 2071 as traced by the H2 1-0 
S(l) 2.12 fim emission ("after Eisloffell l2000h . The driving source of the NGC 2071 outflow 
is believed centered close to RA OS'' 47™ 04?6, Dec +00° 21' 48" (J2000). The box labeled 
Spitzer/IRS bounds the portion of the northern outflow lobe observed here. 
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Fig. 2. Map of the northern lobe of the outflow from NGC 2071 as traced by H2 S(0) through 
S(3) emission. Offsets are relative to the measured peak of the H2 17.035 fim S(l) emission 
centered at RA 05^' 47"^ 08?1, Dec +00° 22' 50.7" (J2000). The peak intensity for each transition 
is provided in Table 1. Positions PI, P2, and P3 are fiducials that mark the positions of the first, 
second, and third strongest pe aks in the S(l) em ission. The position of other known sources 
in the field are also shown (cf. IWalther et a/.lll993|) . 



Page 23 




20 -20 -40 -60 20 -20 -40 -60 

RA Offset (arcsec) RA Offset (arcsec) 



Intensity/ Peak Intensity 



1 I I I I I I I I I f 
0.2 0.4 0.6 0.8 1.0 

Fig. 3. Same as Fig. 2, except for H2 S(4) through S(7) emission. 
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Fig. 4. Same as Fig. 2, except for H2O 725— 616 29.837 jim, 634—505 30.899 fim (blended with 
H2O 854—743 30.871 fim), and 734—625 34.549 /im emission. The H2O 734—625 map is trun- 
cated because Spitzer/IRS was exposed to a strong source while slewing to the NGC 2071 field. 
This caused latent signal levels at the long wavelength end of LH and, thus, we conservatively 
omitted these less reliable observations. 
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Fig. 5. Same as Fig. 2, except for HD R(3) 28.502 /xm and R(4) 23.034 /xm pure rotational 
emission. 
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Fig. 6. Same as Fig. 2, except for SI ^Pi -^Ps 25.249 /xm, Fe II 0*^07/2 - a^Dg/s 25.988 /xm, and 
Sill ^P3/2-^Pi/2 34.815 fim emission. 
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Fig. 7. Top: Contours of H2O 725— 616 29.837 /im emission superposed on the lEisloffell (|2000|) 
H2 1-0 S(l) 2.12 /im map. Middle: Contours of H2 S(l) emission superposed on the map of 
H2O 725— 616 29.837 /im emission. Bottom: Contours of H2 S(l) emission superposed on the 
map of HD R(3) 28.502 /^m emission 
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Fig. 8. Portion of the Spitzer/IRS spectra obtained in a 15" synthesized beam toward H2O Peaks 
1 and 2 (see text). The H2^^0 transition(s) associated with each detected feature are shown. In 
the absence of higher spectral resolving power, we are not able to observationally determine 
the extent to which each transition contributes to the total flux measured for each blended fea- 
ture. In addition to Si II and S III lines in the band, we also note the tentative detection of OH 
2n3/2 19/2+-17/2- 30.277 /xm, ^Hg/s 19/2--17/2+ 30.346 /im, ^Hi/s 17/2+-15/2- 30.657 /im. 
Unidentified lines are denoted with a "U". 
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Fig. 9. Portion of the Spitzer/IRS spectra obtained in a 15" synthesized beam toward H2O Peaks 
1 and 2 (see text). Of particular note is the ground rotational transition of para-H2 at 28.219 /im, 
two HD rotational transitions, R(3) and R(4), at 28.502 and 23.034 /im, respectively, and the 



tentative detection of OH ^Ilg/a 21/2" 



19/2+ 27.393/im, 2n3/2 21/2^ 



19/2- 27.454/im,2n 



1/2 



19/2- - 17/2+ 27.652 /xm, ^Hi/a 19/2+ - 17/2- 27.697 /im, ^Hg/s 23/2+ - 21/2- 25.035 /im, and 
^03/2 23/2- - 21/2+ 25.090 /im. Unidentified lines are denoted with a "U". 
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Fig. 10. Portion of the Spitzer/IRS spectra obtained in a 15" synthesized beam toward H2O 
Peaks 1 and 2 (see text). Of note are the prominent lines of H2 S(l) through (7). Both spectra 
also show broad PAH emission features with peak wavelengths of 6.22, 7.41, 7.59, 7.85, 8.61, 
and 1 1.3 yum. 
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Fig. 11. Root mean square goodness-of-fit of the computed fluxes at each density for H2 
S(0) — S(7) {dashed line-solid squares); HD R(3) and R(4) {dashed line-solid circles); H2O 
Case 1 - the H2O features listed in Table 2, with the exception of the 34.55 /im feature, and as- 
suming an ortho-to-para ratio of 3: 1 {open squares); H2O Case 2 - the pure ortho- or para-H20 
lines at 29.84, 33.13, and 36.21 //m and assuming an ortho-to-para ratio of 3: 1 {solid diamonds); 
H2O Case 3 - the H2O features listed in Table 2, with the exception of the 34.55 /im feature, and 
assuming an ortho-to-para ratio that has been allowed to vary to produce the best fit to the data 
at each density {open circle); H2O Case 4 - the pure ortho- or para-H20 lines at 29.84, 33.13, 
and 36.21 /im and assuming an ortho-to-para ratio that has been allowed to vary to produce the 
best fit to the data at each density {solid triangles). 
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Fig. 12. Total (ortho+para) water abundance, at each assumed density, that produces the best 
overall fit to: (1) the H2O features listed in Table 2, with the exception of the 34.55 /xm feature, 
and assuming an ortho-to-para ratio of 3: 1 {dashed line); (2) the H2O features listed in Table 2, 
with the exception of the 34.55 jim feature, and assuming an ortho-to-para ratio that has been 
allowed to vary to produce the best fit to the data at each density {dot-dashed line); (3) the 
29.84, 33. 13, and 36.21 fim features and assuming an ortho-to-para ratio of 3: 1 {solid line); and, 
(4) the 29.84, 33.13, and 36.21 jim features and assuming an ortho-to-para ratio that has been 
allowed to vary to produce the best fit to the data at each density {dotted line). 
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Fig. 13. Top: Observed spectrum within 15" synthesized beam toward H2O Peak 1. Bottom: 
Model spectra assuming n(Yi2) = 10^ cm^'' and an H2O ortho-to-para ratio of 3:1 (solid line) 
and an H2O ortho-to-para of 1.55:1 (dotted line), which provides the best fit to the pure ortho 
29.84 fim and pure para 33.13 and 36.21 /im lines at this density. 
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Fig. 14. Computed ^^CO line intensities expected for a CO abundance of 10^"^ relative to H2 (see 
text). The dashed line corresponds to the Herschel/FACS spectral line sensitivity of ~ 2x 10^^ 
erg cm^^ s^^ sr (5cr in 1 hour, corresponding to a flux of 5x 10^^^ W m"^ into a 9.7x9.7" 
pixel). 
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Fig. 15. Computed line intensity ratios of select CO high- J rotational transitions with respect 
to the ^^CO J = 15— 14 line measurable with the Herschel Space Observatory. These line ratios 
are particularly sensitive to H2 densities between 10^ and about 3 x 10^ cm^^. 



